A review on recent results on on/off reconstruction of optical states 



G. Brida 1 , M. Genovese 1 , M. Gramegna 1 , A. Meda 1 , S. Olivares 23 , 
M. G. A. Paris 3 2 4 , F. Piacentini 1 , E. Predazzi 3 , P. Traina 1 

1 I.N.RI.M., Strada delle Cacce 91, Torino, Italia 
2 CNISM UdR Milano Universitd, 1-20133 Milano, Italia 
3 Dipartimento di Fisica dell 'Universitd di Milano, 1-20133 Milano, Italia 
4 Institute for Scientific Interchange, 1-10133 Torino, Italia and 
5 INFN and Dipartimento di Fisica Teorica dell 'Universitd di Torino, via P. Giuria 1, 10125 Torino, Italia 

The knowledge of the density matrix of a quantum state plays a fundamental 
role in several applications, from quantum information to Standard Quantum 
Mechanics foundations investigation. 

Recently, a method has been implemented in order to obtain the reconstruction 
of the (diagonal) elements of the density matrix exploiting the information 
achievable with realistic on/off detectors (e.g. silicon avalanche photo-diodes). 
Purpose of this paper is to introduce this on/off method and present the last 
experimental results obtained with it in several different cases. 

© 2008 Optical Society of America 
OCIS codes: 42-50-p, 03.67.Hk, 42.62.Eh 



1. Introduction 

The knowledge of density matrix of a quantum state 
is fundamental for several applications, ranging from 
quantum information^ to the foundations of quantum 
mechanics^ and quantum optics^* 5 * 6 * 7 * 8 * 9 * 1 211U. Many ef- 
forts have been made to find reliable methods to fully or 
partially reconstruct the density matrix especially in the 
photon number basis. In this last case, the reconstruc- 
tion of the diagonal elements, i.e., the photon statistics, 
is of extreme relevance for a first characterization of the 
state. Nevertheless, the choice of a detector with internal 
gain suitable for the measurement is not trivial when the 
flux of the photons to be counted is such that more than 
one photon is detected in the time-window of the mea- 
surement, which is set by the detector pulse-response, 
or by an electronic gate on the detector output, or by 
the duration of the light pulse. In this case, we need 
a congruous linearity in the internal current amplifica- 
tion process: each of the single electrons produced by 
the different photons in the primary step of the detection 
process (either ionization or promotion to a conduction 
band) must experience the same average gain and this 
gain must have sufficiently low spread. The fulfillment 
of both requisites is necessary for the charge integral of 
the output current pulse be proportional to the number 
of detected photons. Photon detectors that can operate 
as photon counters are rather rare. Among these, Photo- 
Multiplier Tubes (PMT's) 13 and hybrid photodetectorsP 
have the drawback of a low quantum efficiency, since the 
detection starts with the emission of an electron from 
the photocathode. Solid state detectors with internal 
gain, in which the nature of the primary detection pro- 
cess ensures higher efficiency, are still under development. 
Highly efficient thermal photon counters have also been 
used, though their operating conditions are still extreme 



(cryogenic conditions) to allow common us e 1 15 * 16 ! The 
advent of quantum tomography provided an alternative 
method to measure photon number distributions^. How- 
ever, the tomography of a state, which has been applied 
to several quantum states^, needs the implementation of 
homodyne detection, which in turn requires the appro- 
priate mode matching of the signal with a suitable local 
oscillator at a beam splitter. Such mode matching is a 
particularly challenging task in the case of pulsed optical 
fields. 

Photodetectors that are usually employed in quantum 
optics such as Avalanche Photo-Diodes (APD's) oper- 
ating in the Geiger mode^U^ seem to be by definition 
useless as photon counters. They are the solid state pho- 
todetectors with the highest quantum efficiency and the 
greatest stability of the internal gain. However, they have 
the obvious drawback that the breakdown current is inde- 
pendent of the number of detected photons, which in turn 
cannot be determined. The outcome of these APD's is 
either "off' (no photons detected) or "on" i.e. a "click", 
indicating the detection of one or more photons. Actu- 
ally, such an outcome can be provided by any photode- 
tector (PMT, hybrid photodetector, cryogenic thermal 
detector) for which the charge contained in dark pulses 
is definitely below that of the output current pulses cor- 
responding to the detection of at least one photon. Note 
that for most high-gain PMT's the anodic pulses corre- 
sponding to no photons detected can be easily discrim- 
inated by a threshold from those corresponding to the 
detection of one or more photons. 

It may thus appear sur prising at fi rst that, among the 
various theoretical studies^ 1 - 9 * 20 * 21 * 22 * that have been ad- 
dressed to achieve the reconstruction of the (diagonal) 
elements of the density matrix exploiting the information 
achievable with realistic dete ctors, a favo urable experi- 
mental test has been found irl 23 | 24 | 25 | 26 | 27 l where a high- 
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fidelity reconstruction of the statistics of mono-partite 
and bi-partite quantum optical states was obtained by 
using on/off detector j 21 * 22 l 

Purpose of this paper is to introduce this on/off 
method and present the last experimental results ob- 
tained with it. 



2. Reconstruction method: diagonal elements of 
monopartite systems 

Let us consider a single-mode quantum optical state. All 
the accessible information on the state is expressed in 
its density matrix p, which in the photon number basis 
reads: 
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In particular, the information regarding photon distri- 
bution of the state is given by the diagonal elements 
p n = p nn of this matrix. 

In this section we are going to show how reconstruction 
of the pnS for a general quantum optical state is possible 
via on/off detectors. 

Assuming that our state is to be revealed by a two- level 
detector, like silicon APDs (avalanche photo-diodes) or a 
photomultiplier operating in Geiger mode, i.e. a detec- 
tor discriminating only between vacuum, "0" , and light 
passage, "1", with quantum efficiency < 77 < 1, the 
measurement process is described by a two- value positive 
operator-valued measure (POVM): 
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where the coefficients matrix V (V 77^ : i ^ j =>■ r]i ^ rjj) 
is a nonsingular Vandermonde matrix of order N + 1. 
Since the p array reconstruction via matrix inversion 
would require a huge and almost impossible to obtain 
number of experimental runs to avoid very large numer- 
ical fluctuations, the need of a faster solution protocol 
leads to consider eq.s ^ and Q as a statistical model 
for the p n coefficients, presenting a linear positive prob- 
lem to be solved by means of the maximum likelihood 
(ML) method 28 . By using the expectation maximization 
(EM) algorithm, after imposing the normalization con- 
straint Y] p n = 1, we arrive at the iterative formula: 
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with: 



p>n = reconstructed p n at the i-th iteration; 



non-click events fraction; 



evaluated non-click probability 



(7) 



with the reconstructed p n 



where n 0u and n n are respectively the registered non- 
click events and the total exper imental runs performed 
with quantum efficiency ?7 / J 20 * 2 H 

The confidence interval on the determination of the 
clement p$ can be estimated in terms of the variance: 
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being F n the Fisher information on the reconstructed el- 
ement, given by the formula: 



The non-detection probability, po — Tr[pIIo(?7)], is thus 
given by: 
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Po(v) = 5^(1 - r,Tp n = J2 Mv)Pn, (3) 

ra=0 n=0 

and the detection probability is pi (77) = 1 — po (77) . Sup- 
posing to have the possibility of varying the quantum 
efficiency 77 of our detector, we can perform K measure- 
ments on our optical state, each one with a different 77, 
building the ensemble: 
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where we define A^ n = A n (j] n ). 

If we restrict ourselves to the approximation p n ~ 
Vri > N and we are able to perform at least K different 
measurements (K > N + 1), for the non-click probability 
ensemble ^ we can write the linear system: 
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where: 
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To check the convergence of this algorithm, we look at 
the error parameter e'*), defined as: 
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stopping the iteration number when the value of goes 
below a certain threshold. 

Finally, the fidelity of our reconstructed photon statis- 
tics with respect to the theoretical one, p n the °\ is given 
by the formula: 
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Several tests on this protocol have been made, all of 
them successfully achieving a good agreement between 
the reconstructed p n and the expected ones. As a first 
example, we consider a weak coherent state generated by 
a cw He-Ne laser emission. In fig [l] is possible to view the 
corresponding reconstructed photon number distribution 
(compared to the theoretical one). A best fit procedure 
shows that the reconstructed photon number distribution 
is compatible with the one expected for a coherent state 
with a mean number of photons On) =5.39 
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Fig. 1. Reconstruction of the photon distribution for the 
weak coherent state (blue dots) compared to simulated 
Poissonian distribution (red dots). The data are com- 
patible with the predictions for a coherent state with a 
mean number of photons (n) — 5.39 

As a second example, we mention the experiment on 
single-photon state s generated by producing PDC her- 
alded photons^ESwhere, a pair of correlated photons (of 
different polarization) was generated by pumping a type- 
II /3-barium-borate (BBO) crystal with a CW argon ion 
laser beam (351 nm) in collinear geometry. After having 
split the photons of the pair by means of a polarizing 
beam splitter, the detection of one of the two by a silicon 
avalanche photodiode detector (SPCM-AQR-15, Perkin 
Elmer) was used to herald the presence of the second pho- 
ton in the other channel (namely, a window of 4.9 ns was 
opened for detection in arm 2 in correspondence to the 
detection of a photon in arm 1); this "heralded" photon 
was then measured by another APD (SPCM-AQR-15, 
Perkin Elmer) preceded by an iris and an interference fil- 
ter (IF) centered at 702 nm (4 nm of FWHM) inserted 
with the purpose of reducing the noise due to the stray 
light. The quantum efficiency of the detection apparatus 
was r\ m ax = 20%: lower quantum efficiencies were simu- 
lated by inserting calibrated neutral filters (NF) on the 
optical path. 

The reconstructed statistics perfectly agreed with the 
expected ones: together with a dominant single photon 
component, a vacuum and double pair components were 
measured, respectively of (2,7±0,2)% and ((1,9±0,2)% 
(in agreement with the theoretical values). Finally, we 
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would like to acknowledge a recent measurement where 
the method was applied to stimulated PDC 29 . In this 
case we used a new algorithm where a constraint (on 
photon number) was introduced (through Lagrange mul- 
tipliers), showing as this allows a very good reconstruc- 
tion in critical cases. 

3. Extension to the bipartite case: the theory 

Since for many applications multipartite states are 
needed, our aim here is to extend our model of the previ- 
ous section to this case as well. In particular, for the bi- 
partite case (easily generalizable to the multipartite one) , 
we will consider two propagation modes of our optical 
field (i.e. the two outputs of a beam splitter) and study 
both of them at the same time, tryin g to extend the pre- 
vious formulas to this new frameworlPS 
To begin with, instead of the probabilities (J3j) we will 
have the joint click/no-click probabilities: 



(12) 



where: 

Qnk = (nk\g\nk) (\nk) = \n) ® \h)\ (13) 

is the joint photon distribution of the state and r\ the 
quantum efficiency of the hypothetical detectors consid- 
ered. 

The above equations provide a relation between the 
statistics of the clicks of our detectors and the actual 
statistics of photons. Again, let us restrict our consider- 
ations to the case g nk ~0Vn,fc>iV(i.e. we "confine" 
our bipartite state into a (AT+1) x (AT+1) Hilbert space): 
if we can properly change the quantum efficiency of our 
system in such a way that K different measurements can 
be performed (with K different values /i — 1, ...,K, 
ranging from rji — rj m i n to a maximum value r/x — ?7max)> 
the whole amount of on/off detection statistics collected 
can give enough information to reconstruct the joint pho- 
ton distribution of the bipartite state. 

In details, by defining the vectors: 

g=(p&,... i p%,&,...,p$,pr Q ,...,p%) (w) 

q = (Qoo, Qoi, Qw, — ) (15) 
Qnk^q P (p = l + k + n(l + N)), (16) 

we can summarize the on/off statistics with the compact 
formula: 
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where we have introduced the matrix B with entries: 
A^A^k if (i = l,..,K 

,,, \ Apnil — A^k) if n = K+l,..,2K (18) 

(l-A^A^ if M = 2X+l,..,3ir 

where k = (p- l)mod(l + iV) and n = (p-l-k)/(l+N). 
In our conditions, eq. ( 17 1 represents a finite statistical 



linear model for the positive unknown q p . Like in the 
monopartite case, by means of the maximum-likelihood 
(ML) principle we achieve a well approximated solution 
of this LINPOS problem with the iterative algorithm 28 : 
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In eq. ( 19 1 we have: 

(7p = p-th element of the reconstructed photon 
statistics at the i-th iteration; 
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^Qp^~\ — evaluated non-click probability 
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from the reconstructed q v 
/i = K + 1,..,2K 
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being noi?7 , nio^ > ^oor^ the number single and double 
non-click events observed on the whole amount n p of ex- 
perimental runs made with r; = rj^. 

To evaluate the confidence interval or. 
tion of the element qffl we still use eq. ( i 
Fisher inforr 
rewritten as: 
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The analogous of the total error v- % > written in eq. ( 10 ) 
is given by: 
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which measures the distance of the reconstructed non- 
click probabilities from the measured ones: the algorithm 



is stopped when reaches its minimum, or goes below 
a certain threshold value. 



Finally, the fidelity equation (111 for the reconstructed 



joint photon number distribution takes the form: 
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giving us the chance to compare the obtained q p with 

the expected ones {qp heo ^). In the end, here we have pre- 
sented our new photon statistics reconstruction method 
(based on on/off detection) for a generic bipartite optical 
field: the next section will be dedicated to the experimen- 
tal tests we made in order to check its reliability. 

4. Extension to the bipartite case: experimental 
tests 

Let us now report the experimental results obtained by 
applying the reconstruction method presented above to 
two different situations, the single photon state passing 
through a beam splitter (BS) and the split PDC single 
branch: as we will see soon, our algorithm managed to 
obtain a good reconstr uction of the joint photon distri- 
bution in both the case d 26 * 2 ^ . 

A. Setup A: Type-II PDC single photon source 



In our first experimental setup (fig. |2pl a 0.2 W, 398 
nm pulsed (with 200 fs pulses and 70 MHz repetition 
rate) laser pump was generated by second harmonic of 
a Ti:Sapphire laser at 796 nm; it was then injected into 
a 5 x 5 x 1 mm type-II BBO crystal, generating entan- 
gled photon pairs by Parametric Down-Conversion. The 
detection of a photon on one of two correlated branches 
of degenerate PDC emission was then used as trigger to 
herald the presence of the correlated photon in the other 
direction. 

Thus, the idler photon was addressed to an optical filter 
(a narrow band interference filter or a red glass, as will 
be specified later), collected by a lens and then sent to 
a silicon APD (APD1). Instead, the corresponding sig- 
nal was properly filtered (with the same filter as the idler 
photon) and then injected into a beam splitter (BS), sep- 
arating in two the optical path of the photon and thus 
generating a bipartite state that is the superposition of 
a single photon (|1)) with the vacuum state (|0)): 



l^Bs) = \^|oi) + \/r^|io) 



(24) 



where r is the BS transmittance. 

The BS was followed on both output arms by a collec- 
tion/detection apparatus, respectively called APD2 and 
APD3, of the same type as APD1 (all the detectors were 
Perkin Elmer SPCM-AQR-15 silicon APDs). 

The proper set of quantum efficiencies was obtained 
by inserting before the BS several neutral filters (NF) 
of different transmittance, calibrated by measuring the 
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Fig. 2. schematics of the experimental setup featuring the 
Type-II PDC heralded photon source. The idler photon 
is addressed to an IF (RG) filter, collected and sent to 
APD1, opening a coincidence window in the TAC mod- 
ules; the signal goes through the NF and the IF (RG) 
filters, and then is split by the BS, whose outputs are 
collected and sent to APD2 and APD3 to close the co- 
incidence windows opened. The output of the two TACs 
is also sent to an AND logical gate whose outputs gives 
the number of double coincidences. 



ratio between the counting rates on D2 and D3 with the 
filter inserted and without it. 

In correspondence of the detection of a photon in arm 
1, a coincidence window was opened on both detectors on 
arm 2: this was obtained by sending the output of Dl as 
Start to two Time-to- Amplitude Converters (TAC) that 
received the detector signal of D2 and D3 as Stop. The 
20 ns coincidence window was set such to avoid spurious 
coincidences with PDC photons belonging to the follow- 
ing pulse (we remind that the repetition rate of the laser 
was 70 MHz). 

The TAC Out outputs were then addressed both to 
the computer and to an AND logical gate in order to 
reveal coincidences between them; its output was also 
collected via computer, together with one TACs Valid 
Start (giving us the total number of coincidence windows 
opened). 

These four data sets allowed us to evaluate the 
click/non-click frequencies /oq, fox, fxo, fix, needed for 
the reconstruction of the photon statistics of the bipartite 
state with our recursive method (except for fix)- The 
background is estimated and subtracted by measuring 
the TACs and AND outputs out of the window triggered 
by APD1 detection; then, the maximum quantum effi- 
ciency r\ m ax was calculated as the ratio between the sum 
of coincidences in APD2 and APD3 and th e cou nting rate 
on APD1 without the insertion of any NF 31 | 32 [ 

With the purpose to verify the method in different 
cases we considered four alternatives given by the com- 
bination of a balanced (50% - 50%) or unbalanced (40% 



- 60%) BS with two optical filters sets, either large band 
red glass filters (RG) with cut-off wave length at 750 nm, 
or interference filters (IF) with peak wave length at 796 
nm and a 10 nm FWHM. 

The first test was made with the 50%-50% BS and the 
interference filters: we collected data for K = 33 different 
quantum efficiencies. Elaborating these data with our re- 
construction algorithm within a 3 x 3 Hilbert space choice 
(JV = 2) lead to the reconstructed joint photon distribu- 
tion shown in fig. [3j here we can appreciate how the only 
relevant entries are qqx and gxo (single photon transmit- 
ted or reflected by the beam splitter), well reproducing 
the inserted BS ratio. 




Fig. 3. reconstructed g n k entries of the joint photon dis- 
tribution of our \ipBs) for the setup with 50%-50% BS 
and 10 nm FWHM IF: within the uncertainty intervals 
we can see poi — Sxo> m agreement with the BS r value. 

There is also a low but non-null vacuum component 
qoq: this can be due to a non perfect background eval- 
uation and subtraction in the experimental data, or to 
some light absorption in the system (maybe caused by 
the BS cube). 

The second test was performed with the same BS but 
with different optical filters, namely red glass (RG) fil- 
ters with X cut _ ff = 750 nm: the maximum quantum 
efficiency obtained this time was rj max = 8, 8%, and the 
acquisition was again repeated K = 33 times. Similar 
results were obtained with red glass filters. 

Then we replaced the 50%-50% BS cube with an un- 
balanced 60%-40% BS plate, maintaining the RG large 
band filters: with this setup we had r\ max — 12,3%, and 
performed a K — 41 steps data collection. The difference 
between the transmitted and reflected branch of the BS 
was evident and in agreement with the known beam split- 
ter ratio, also we did not have any non-null g n k compo- 
nent except for gxo and goi, thus the reconstructed state 
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is definitely the one of eq. (24 1 with t = 0,4. Finally, for 



the last test we replaced the RG filters with the previous 
IF ones obtaining a similar result. 

B. Setup B: Type-I PDC thermal source 

As a second example we considered a single branch of 
PDC emission without triggering, corresponding to a 
multi-thermal statistics^. 

In a setup similar to the previous onep^, we generated 
PDC light by means of a 5 x 5 x 5 mm Type-I BBO 
crystal pumped by a Q-switched (triplicated to 355 nm) 
Nd:Yag laser with 5 ns pulses, power up to 200 mJ per 
pulse and 10 Hz repetition rate. Because of the very high 
power of the pump beam, a state with a large number 
of photon was generated each pulse. We have therefore 
attenuated (by using 1 nm FWHM IF and neutral filters) 
the multithcrmal state before detection. 

Again, the different quantum efficiencies were obtained 
by inserting (before the BS) Schott neutral filters, whose 
calibration this time was obtained by measuring the 
power of a diode laser before and after them, with the 
calibration laser injected in the same point as the PDC 
in order to minimize the effect due to eventual non homo- 
geneous NF filters. The coincidence scheme was realized 
sending two Q-switch triggered pulses to two TAC mod- 
ules as Start inputs, and the APDs outputs as Stop; then, 
having set properly the 20 ns coincidence window, we ad- 
dressed the two TAC Out outputs to the AND logic port, 
and the Valid Stops to the counting modules (together 
with one TAC's Valid Start and the AND output). The 
background was estimated and subtracted in the same 
way as in the previous experiment. 

The maximum quantum efficiency f} max was in this 
second case evaluated by multiplying the APDs nomi- 
nal quantum efficiencies, the IF peak transmittance and 
the fiber couplers effective efficiency, measured with the 
diode laser: it turned out to be r\ max — 25,0%. 

The expected on/off joint statistics for this optical 
state is: 



In summary, we can say that this experiment gave an- 
other confirm of the reliability of our algorithm, even 
when dealing with a larger Hilbert space. 

5. Reconstruction of the non-diagonal elements 

The previously described method can deliver complete in- 
formation on a quantum optic state when the density ma- 
trix of such a state is diagonal, and all the non-diagonal 
elements are equal to zero. Obviously, a more general 
procedure is required which takes into account also the 
states in which not all the non diagonal elements vanish. 
In order to reconstruct also the non-diagonal elements, 
and thus the entire density operator, it will be shown that 
additional phase information is needed. Before entering 
into the details of our implementation, in order to intro- 
duce the theoretical method^, let us just suppose for the 
moment that the state to be reconstructed is mixed with 
a local oscillator mode by an unbalanced beam-splitter 
(with high transmittance and low reflectance). It is easy 
to verify that, in this case, the transmitted mode is equiv- 
alent to the signal mode shifted by a displacement oper- 
ator D(a) — exp(aa^ — a* a), where a (a') is the photon 
destruction (creation) operator associated to the signal 
mode and a — \a\ e %H> is the local oscillator field am- 
plitude. Measuring the photon statistics of the output 
state is equivalent to measuring the displaced Fock-state 
probability distribution of the signal: 
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where \n, a) = D(a)\n) are the displaced Fock states and 
p is the signal mode density operator. Note that p n is 
now a function of the displacement a also. 

For any physical state the density matrix elements 
(fe|p|m) decrease with increasing k and m, therefore 
p n (a) can be approximated to any desired degree of ac- 
curacy by setting those elements to zero for k,m > n , 
where no is & conveniently large integer number. 

It follows that we can expand eq. [27] in the form: 
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p\ = {M[M + V {1 r)^.]" 1 }" - pg 



where N ave is the average number of photons, M the 
number of propagation modes and r the BS transmit- 
tance. 

The reconstructed joint photon statistics (now upon a 
17 x 17 truncated Hilbert space) was compared with the 
multithcrmal distribution: 



(n + m + M-l )\ (1 + %^) 
n\m\{M-l)\ 
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The calculated fidelity, with the formula (11), for i > 
2000 was always greater than 99%. 
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Expressing the displaced Fock states \n, a) in the ordi- 
nary Fock basis, one obtains^! 



p n (a) = e ! Q ' 2 n! ^ v / k\m\(k\p\m)x 

k,m— 

(_l)j+J| a F+'+2(H-l) e i( m -k) f 



3=0 1=0 



j\(n - j)\{k - j)W(n - l)\(m - /)! 



where j = min{n, k} and I = min{n, m}, so that for any 
value of | a |, p n {cx,) can be regarded as a function of tp 
and expanded in a Fourier series, the general component 
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being: 
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with: 

G« ro (|a|) = 
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(31) 



where j = min{n,m + s}. We notice that p^dal) is 
related to the density matrix element whose row and 
column indices differ by s. If the photon number dis- 
tribution Pn(ct) is measured for n — 0,1,..., N, with 
N > n ,eq. ( |30| for each s represents a system of (N+ 1) 

linear equations linking (N + 1) measured quantities pn^ 
to (no + 1 — s) parameters (the unknown density matrix 
elements). This system is clearly overdetermined, so it 
can be inverted using the least squares method in order 
to obtain the density matrix elements from the measured 
probabilities. 

The reconstructed off-diagonal density matrix ele- 
ments can thus be obtained as: 

N 

(m + s\p rec \m) = £F(t(H)^)(H) (32) 



?i=0 



with: 



41(H) = {[Gll(|a|)] T G( s i L(|a|)}- 1 [Gil(|a|)] T 

(33) 

The F matrix satisfies the condition: 



iV 



E F S»(N) G "t(W) = <w 



(34) 



n=0 



for m, ml = 0, 1, . . . , n — s, so that from the exact prob- 
abilities the correct density matrix elements are found 
(p rec = p). Furthermore, the least squares method en- 
sures that the *pffl calculated from p rec according to eq. 
5asured 

N 

Eoi s) 



(30 1 best fits the measured quantities such that: 

N 



P 



(35) 



is minimized. 



In conclusion, combining eq.s (32) and (30 1, we find 



out the formula for the direct sampling of the density 
matrix from the measured photon number distribution 
of the displaced stated 



(m+ s\p rec \m) 



N 

^E/ Fi%{\a\)e is *p n { a )dv (36) 

n=0 



It is worth mentioning that this method, given that the 
p n (a) are known, requires only the value of tp to be var- 
ied. 

Let us now extend our method for the case of non- 
perfect detection, in which the detection efficiency rj is 
less than unity. The effective measured photocount dis- 
tribution Pk(a) is related to the photon number distri- 
bution p n {ot) as follows: 



n=0 



where Mk in (rf) is 
M k ,n{ri) = 



(37) 



( n k )r] k {l -r]) n - k k<n 
k > n 



(38) 



since rj k (l — i]) n ~ k is the probability to detect k among 
n incoming photons and (?) is the number of possible 
k- photons subsets. So, the equivalent of eq.(30l for the 
measured quantities Pi s ^(|a|) yields: 

n — s 

4 S) (M)= E G^l(\a\,r]){m + s\p\m), (39) 

m=0 

where the new matrices: 



Gi%(\a\, V ) = Y,MnAv)G[ s > m (\*\) 



(40) 



fc=0 



obtained from the Gj. s ^j(|a|) defined in eq. (31), can be 
inverted in the same manner described above to ob tain 
some matrices Fm^ n (\a\ ,77) to be used, as in eq. (36 1, to 
reconstruct the density matrix. 

6. Experiment: full density matrix reconstruc- 
tion for a coherent state 

The key element for the reconstruction of non-diagonal 
elements in the density matrix of the measured state is 
a sort of asymmetric homodyne detection in which the 
signal mode and a local oscillator are mixed by an unbal- 
anced beam splitter. In our scheme, actually, the local 
oscillator is not an external source but it is obtained by 
the same source as the signal divided in two by a 50/50 
beam splitter. 

To be more precise, the source beam enters in a Mach- 
Zchnder interferometer in which the part reflected by the 
first beam splitter is taken as the signal, as it travels 
along the "long" path inside the interferometer, while 
the transmitted portion is regarded as the local oscilla- 
tor. The two modes are then mixed by the second beam 
splitter, which is effectively used to perform a sort of un- 
balanced homodyning. It is possible to change the length 
of the optical path in the long arm of the inteferometer, 
thus varying the phase or the signal mode with respect 
to the reference mode. 



The expression of the mean number of photons at the 
output port of the interferometer it is easy to be calcu- 
lated and yields: 



(n ) = Tr[a\a p] = {n, 



1 — cos(ip) 



(41) 



where p and hi are, respectively, the density and the 
photon number operators for the input state. We re- 
mind here that the non-click probability po, which is the 
quantity to be actually measured, for a coherent state of 
amplitude v, reads: 



fa) = I>-^rj>n = I>-'7) B - 



2n g->)M 2 



(42) 



n=0 



which obviously is a function of the overall quantum ef- 
ficiency i] of the detection system. In our case, r>n_ is 



41 



also a function of tp since \v\ can be derived from eq 
(H 2 = (n)). 

The theoretical method exposed in the previous section 
leads to the possibility of a direct experimental sampling 
of the density operator of the initi al stat e by performing 
photon distribution reconstructions^^ f t ne state exit- 
ing the interferometer for a suitable number of tp values. 

We point out that our scheme allows to estimate the 
error in real time, that it is suitable both for pure and 
mixed states and that statistic information is obtained by 
sampling a discrete matrix rather than measuring a con- 
tinuous distribution in phase space as homodyne quan- 
tum tomography schemes. 

A. Setup 

In our set-up (see figQthe signal state is the output of 
a He-Ne laser (A = 632.8 nm), which is lowered to sin- 
gle photon regime by neutral filters. The spatial profile 
of the signal is purified from non-gaussian components 
by a spatial filter realized by two converging lenses and 
a 100/im diameter wide pinhole. An iris just after the 
pinhole ensures the selection of a single Gaussian spatial 
mode. The laser cavity is also preserved by backreflec- 
tions, which may cause instability, by means of an optical 
isolator consisting in a Faraday rotator between two po- 
larizers. The second polarizer (let us name it B) angle 
is shifted of 45 degrees with respect to the first (A), so 
that the light transmitted by the latter, whose polariza- 
tion is rotated by 45 degrees by the Faraday rotator, is 
all transmitted also by the second. Since the polariza- 
tion rotation by Faraday effect is in the same direction 
regardless of the laser propagation direction, any backre- 
flected light passing through B, would suffer another 45 
degrees rotation before reaching a polarizer and would 
thus be stopped, being orthogonal to the polarizer angle. 

After a beam-splitter, part of the beam is addressed 
to a control detector in order to monitor the laser ampli- 
tude fluctuations, while the remaining part is sent to the 
interferometer, its main structure consisting in a single 
invar block custom designed and developed at INRIM. 




0<%<66% 



Fig. 4. Setup for the reconstruction of the density ma- 
trix for a coherent state. The emission of a He-Ne laser 
(X = 632.8 nm) is lowered to single photon regime by 
neutral filters. A spatial filter realized by two converg- 
ing lenses and a 100/im diameter-wide pinhole purifies 
the shape of the signal and allows to select a single spa- 
tial mode. A beam-splitter reflects part of the beam to 
a control detector used to monitor the laser amplitude 
fluctuations, while the remaining part is sent to the in- 
terferometer. The phase between the "short" and "long" 
paths in the interferometer can be changed by driving the 
position of the reflecting prism by means of a PI piezo- 
movement system. A set of variable neutral filters allows 
to collect photons for different values of the quantum effi- 
ciency. The detectors used are Perkin-Elmer Single Pho- 
ton Avalanche Photodiode(SPCM-AQR) gated by a 20 ns 
wide time window with (repetition rate — 200 kHz). A 
single run consists of 5 repetitions of 4 seconds acquisi- 
tions and events are recorded by a NI-6602 PCI counting 
module. 



A PI piczo-movemcnt system allows to change the phase 
between the "short" and "long" paths by driving the po- 
sition of the reflecting prism with nanometric resolution 
and high stability. 

For each position of the prism, the "no-click" probabil- 
ities are collected for different sets of neutral filters, and 
thus, for different quantum efficiencies. 

The detector, a Perkin-Elmer Single Photon Avalanche 
Photodiode (SPCM-AQR), is gated by a 20 ns wide time 
window with a repetition rate of 200 kHz. In order to 
obtain a reasonable statistics, a single run consists of 5 
repetitions of 4 second acquisitions. Events are recorded 
by a NI-6602 PCI counting module. In this case, since all 
the attenuations in front of the detectors can be included 
in the generation of the state, the overall maximum quan- 
tum efficiency is assumed to be 66.5%, as the nominal 
efficiency declared by the manufacturer data-sheet of the 
photodetectors. 

B. Experimental results 

In our preliminary experiment!^ we consider M v = 11 
number of different phases and reconstruct the corre- 
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sponding photon distributions using the method exposed 
in the previous sections. Due to the small number of 
phases considered, only the diagonal and the first off- 
diagonal elements of the density matrix can be recon- 
structed. The results are shown in Fig. [5] where we 
plot the reconstructed (jo£, e £j, the expected (p^J den- 
sity matrix elements and the absolute difference S = 
\p r n,rn ~ Pn%i\- The reconstructed p rec corresponds to a 
coherent state |a) with energy \a\ 2 = 3.06. Considering 
the small amount of data, the agreement between p rec 
and p exp is good. Another measurement has been per- 



H rec s true g delta 




Fig. 5. In order, from left to right, the three plots show the 
reconstructed density matrix (rec), plfm, for the coherent 
state, the expected matrix (true), p^^, and the absolute 
difference S = |p^ 6 f„ — Prfm\ (delta) between the two. In 
particular, the last plot show good agreement between the 
reconstruction results and the theoretical values in view of 
the small number of phases considered for reconstruction. 

formed with the same setting but without unbalancement 
in the second beam splitter, and the collected data po(^) 
for several values of tp are being evaluated right now. 
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